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Abstract. Given an ordinary differential field K of characteristic zero, it is 
known that if y and 1/y satisfy linear differential equations with coefficients 
in K, then y 1 /y is algebraic over K. We present a new short proof of this 
fact using Grobner basis techniques and give a direct method for finding a 
polynomial over K that y' /y satisfies. Moreover, we provide explicit degree 
bounds and extend the result to fields with positive characteristic. Finally, we 
give an application of our method to a class of nonlinear differential equations. 



1. Introduction 

Using Grobner basis techniques, we provide new constructive proofs of two the- 
orems of Harris and Sibuya jHlE] (see also, H3E] and [HI Problem 6.60]) that give 
degree bounds and allow for several generalizations. Their results were the follow- 
ing. 

Proposition 1.1. Let Ni,N 2 > 1 be positive integers and let K be a differential 
field of characteristic 0. Let F be a (differential) field extension of K and suppose 
that L\{Y) and L 2 {Y) are nonzero homogeneous linear differential polynomials (of 
orders N± and N 2 respectively) with coefficients in K . Further, suppose that one 
of the following holds: 

(1) yGFhas L 1 (y)=L 2 (l/y) = 0. 

(2) N 2 < q e Z+, and y e F has L^y) = L 2 (y q ) = 0. 

Then, Dy/y is algebraic over K . 

In this article, we prove the following more refined result. 

Theorem 1.2. Let N\,N 2 > 1 be positive integers and let K be a differential 
field of characteristic p. Let F be a (differential) field extension of K. Suppose 
that L\(Y) and L 2 (Y) are nonzero homogeneous linear differential polynomials (of 
orders N x and N 2 respectively) with coefficients in K . Further, suppose that one 
of the following holds: 

(1) p > JVi + N 2 - 2 and y e F has L x {y) = L 2 (l/y) = 0. 

(2) p > {q + l)(iVi - 1) - 1, N 2 < q G Z+, and y E F has L x {y) = L 2 {y«) = 0. 
Then, D^y/y is algebraic over K for all j > 1. Moreover, the degree of the poly- 
nomial for Diy/y (j — 1, . . . , Ni — 1) in (1) is at most (^~jyr^ and for (2), at 



1991 Mathematics Subject Classification. Primary 34M15, 13P10; Secondary 34A26. 
Key words and phrases, logarithmic derivative, linear differential equation, differential field, 
Grobner basis. 

This work is supported under a National Science Foundation Graduate Research Fellowship. 



2 



CHRISTOPHER J. HILLAR 



most ^^m^^Y^\ ■ Additionally, if p = 0, then the first condition in (1) and in 
(2) may be dropped. 

This theorem is proved in section 2. In section 3, we give tighter degree bounds 
under the assumption that K has characteristic zero (or more generally, if K is an 
infinite perfect field), and finally in section 4, we describe how our technique applies 
to certain nonlinear differential equations. 

Recall that a polynomial / <E K[x] is called separable if all of its roots are distinct, 
and a field K is called perfect if every irreducible polynomial in K[x] is separable. 
Examples of perfect fields include finite fields, fields of characteristic zero, and, of 
course, algebraically closed fields. It is interesting to note that there is a converse 
to Theorem II .21 for perfect fields K. Specifically, we have the following result. 

Proposition 1.3. Let K be a perfect field. Ify'/y is algebraic over K, then both 
y and 1/y satisfy linear differential equations over K . 

Proof. Suppose that K is perfect and u — y'/y is algebraic over K. Let f(x) = 
x m + a m -ix m + • • • + do 6 K[x] be the monic, irreducible polynomial for u over 
K. Since K is perfect, it follows from basic field theory that gcd(/, ^j) = 1. In 
particular, 7^ 0. Consider now, 

(Tn— 1 \ m— 1 

mu™" 1 + i^u 1 - 1 + ^ a'y. 
i=l / i=0 

Since ^ = mx m ~ x + Y1T=^ iaiX l ~ x is not the zero polynomial, it follows that 
mu™ 1 " 1 + Y^i^ 1 io-iU 1 " 1 7^ (by irreducibility of /). Hence, u' S K(u) and the 
same holds for higher derivatives. 

Next, notice that (1/y) 1 = —y'/y 2 = —u/y and an easy induction gives us that 
(1/y)^ — pk{u,u' ,u" , . . )/y, in which pu is a polynomial (over K) in u and its 
derivatives (set po = 1). By above, the polynomials pk(u,u', . . .) lie in the field 
K(u). This implies that they satisfy some (non-trivial) linear dependence relation, 

JV 

h k p k = 0, 

k=0 

in which h k € K. Therefore, 

N N 
k=Q k=0 

as desired. Performing a similar examination on the derivatives of y' = uy produces 
a linear differential equation for y over K, completing the proof. □ 

As an application of our main theorem, take F to be the field of complex mero- 
morphic functions on C and K = Q. Then, the only y such that both y and 1/y 
satisfy linear differential equations over K are the functions, y = ce ux , in which u 
is an algebraic number of degree at most m\n{Nx, N 2 } and c S C \ {0}. 

This simple example shows that it is possible to produce a minimum degree of 
min{iVi, ^2} for y'/y] however, it is still an open question of whether we can achieve 
a minimum degree close to the bound given in Theorem 11.21 (or for characteristic 
zero, Theorem 13 . II below) . 
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Theorem II . 21 can also be used to show that elements in a differential field F do 
not satisfy linear differential equations over a subfield K, as the following example 
demonstrates. 

Example 1.4. Problem 6.59] Let K = C{x) and F = C((x)). Then, sec(ir) does 
not satisfy a linear differential equation over K. Indeed, cos(x)'/cos(x) = — tan(x) 
is not algebraic over C(a;). 

We begin with a lemma that is fundamental for the solution of this problem, 
although its generality should be useful in many other contexts. Let if be a field 
and let R = K[x±, . . . ,x n ] be a polynomial ring over K. Give R a grading by 
assigning to each Xi, a number w(xi) = Wi G N, so that 



Then, we have the following more general version of a theorem of Sperber [Jj. 
Its proof can be found in Lemma 2.2.2]. 

Lemma 1.5. Let I be the ideal of R generated by a collection of polynomials and 
let J be the (homogeneous) ideal generated by the leading forms (with respect to 
the grading above) of these polynomials. Then, dim(J) > dim(7). 

From this we obtain the immediate corollary that we need. 

Corollary 1.6. If J is zero-dimensional, then so is I. 

2. The Main Theorems 

Before embarking on proofs of the theorems mentioned in the introduction, we 
present a simple example to illustrate our technique. Let y%,y2,--- be variables. 
We will view yj = D J y/y as solutions to a system of polynomial equations over 
^[{l/j }.&.]• F° r example, consider the system (Ni = 3, N 2 = 2): 



in which 02, ai, ao, 6i, &o £ K. Dividing the first equation by y and expanding the 
second one gives us the more suggestive equations: 



Also, differentiating the original equation for X/y and expanding, we have that 



Thus, we may view (y' /y,y" /y,y"' jy) — (2/1,^/2,^3) as a solution to a system of 
three polynomial equations in three unknowns. 

Let w{yi) = i define a grading of K[y±, j/2, 2/3], and notice that the system of 
leading forms, {2/3 = 0, 2y\ — j/2 = 0, — 6y^ + 62/12/2 — 2/3 = 0}, has only the trivial 
solution (2/1,2/2,2/3) = (0,0,0). In light of Corollary 11.61 it follows that the equa- 
tions above define a zero-dimensional variety, establishing Theorem ll.2l (1) for this 
example. 

In general, we will construct a system of N± — 1 equations in A^i — 1 unknowns 
satisfied by the yt. These equations will define a zero-dimensional variety, and thus, 




y'" + a 2 y" + aiy' + a Q y = 0, 
(l/l/)" + 6i(Vl/)' + 6o(Vl/) = 



2/3 + a 2 2/2 + aiyi + a = 0, 
(2y 2 1 - y 2 ) - bxyi +b = 0. 



62/? + 62/12/2 - 2/3 + M22/i 2 " 2/2) - 2/1 (&i + b ) + b' = 0. 
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standard elimination techniques (see, for instance, 1 ) give us a direct method of 
computing, for each i, a nonzero polynomial (over K) satisfied by yi. 

Let us first examine what happens when we compute /„ = D n (l/y). Notice that 

fa = Vy 
fi = -y~ 2 Dy = -yi/y 
h = 2y- 3 (Dy) 2 - y~ 2 D 2 y = 2yf/y - y 2 /y 

fa = —62/1/2/ + 6ym/y - 2/3/2/- 

In general, these functions /„ can be expressed in the form f n — (l/y)p n (yi, . . . , y n ) 
for polynomials p n € Z[yi, . . . , y n ]. More remarkably, with respect to the grading 
w {yi) = h these p n are homogeneous of degree n. These facts are easily deduced 
from the following lemma. 



Lemma 2.1. Let m 6 Z + . Then, 

m — 1 



Pm _ \ P 
m} ^— ' (in 



'm—j yj yn 



ml ^— J [m — j)ljl ml 



Proof. Consider the following well-known identity (Leibniz' rule), 

Q (Wh) (D m ->g) = D m (hg). 



Setting h = y and g = l/y, it follows that 

j\ (m-i)! 

Multiplying the numerator and denominator by y and rewriting this expression 
gives us 

Prn _ Pm~l 2/1 Pm-2 2/2 _ _ Pi 2/m-l _ 2/ro 

ml ~ (m-l)!l! (m-2)!2! "' l!(m-l)! ml' 

□ 

We are now ready to prove Theorem 11.21 ( 1 ) . 

Proof of Theorem \l.H\ (1). With Ni,N2 as in Theorem 11.21 we suppose Ni = n, 
N2 = m. Dividing through by y in the first differential equation for y gives us 

(2.1) 2M = -fln-ij/n-i ai2/i - ao a>i 6 # 

while multiplying the second one for l/y by y produces the equation 

Pm + b m -iPm-i -\ h& = i), e K. 

Differentiating the original linear differential equation for y, k times, we will arrive 
at linear equations y n +k = Lk{y\, ■ ■ . , JM-i) in terms (over K) of y\, . . . , y n -i like 
(|2.1() above (by repeated substitution of the previous linear equations). If we also 
differentiate the equation for l/y, k times, we will produce another equation for 
the variables y^. More formally, we have that 

D m+k (l/y) + D\b m - l D m - 1 {l/y)) + ■■■ + D k (b /y) = 
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produces the equation (by Leibniz' rule) 

tn— 1 k / / \ 

D m+k (l/y) + ££ . (i^M (X> fc - J+< (1/tf)) = 0. 

i=0 j=0 ^ ' 

So finally (after multiplying through by y), it follows that 

(2.2) P m+fc := +££(.) (-D^i) PJfc-i+i = 0. 

i=0 j"=0 

It is clear that the leading homogeneous forms of the P m +k (with respect to the 
grading above) are p m +k- Consider now the ring homomorphism <j> : i^[{?/j}i=i] ~* 
K[yi, . . . , y n -i] defined by sending yj i— > for j > n and j/j i— * yj for j < n. Let 
P m +k denote the polynomials produced by substituting the linear forms Li for the 
variables y n +i (i = 0, 1, . . .) into the polynomials, P m +k- The leading homogeneous 
forms of the P m +k will just be p m +k '■= 4>{Pm+k) because we are substituting linear 
polynomials with strictly smaller degree (corresponding to the grading) . In light of 
Corollary 1 1.61 it suffices to show that the n — 1 equations (in the n — 1 variables), 

(2.3) p m = 0, p m +l = 0, . . . , pm+n-2 = 0, 

are only satisfied by the point (0, . . . , 0) to prove the claim. 

As for the degree bounds in the theorem, simply notice that each polynomial, 
p m +k, has total degree m + k (coming from the y™ +k term), so that Bezout's 
theorem gives us the bound of / +n 7w ■ 

G (m— 1): 

Suppose that (yi, . . . , y n -i) ^ (0, ... ,0) is a zero of the system in l|2.3ll : we will 
derive a contradiction. Let r 6 {1, . . . , n — 1} be the largest integer such that 
y r 0, and choose t S {0, . . . , m — 1} maximal such that p m -t = 0, p m ~t+i = 0, 
. . ., p m = 0. If t = m — 1, then p\ = —y\ = 0, and so the recurrence in Lemma 12. II 
and (2.3) give us that yt = for i 6 {1, . . . , n — 1}, a contradiction. Thus, t < m — 2. 
Using Lemma 2.1 with <f> (and the maximality of r), examine the following identity: 

Pm-t+r-l _ Pm-t+r-2 Hi Pm-t Ur-1 Pm-t-l Vr 



(m-t + r-1)! (m-t + r -2)\1\ (m-t)!(r-l)! (m-i-l)!r! 

From 1)2. 3|) and the property of t above, it follows that ^_~t_~i)! ft = ^" Thus, 
y r = or p m -(t+i) = 0; the first possibility contradicts y r ^ 0, while the second 
contradicts maximality of t. 

This proves that the equations (|2.3|) define a zero-dimensional variety. Finally, 
we remark that our methods still hold if the characteristic of K is larger than the 
denominator (m + n — 2) above. This completes the proof of the theorem. □ 

The proof for Theorem 11.21 (2) is similar to the one above, however, the recur- 
rences as in Lemma 12.11 are somewhat more complicated. Let n S N, q € Z + 
and examine f n>q = D n (y q ). It turns out that f n>q = y q p n ,q(yi, ■ ■ ■ ,Vn) m which 
Pn,q € . . . , y n ] is homogeneous of degree n (with respect to the grading 

w (Ui) — 0- This follows in a similar manner as before from the following lemma. 

Lemma 2.2. Let p n ,\ = y n for neN (j/o = !)• Then, for all m £ N, g > 1, 

Pm,, = y m + £ ( ■ WjPm-^-l- 
3=0 
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Proof. Use Leibniz' rule as in Lemma 12 . II with h = y q 1 and g = y. □ 

The next lemma will be used in the proof of Theorem 1 1.21 (2). and it follows from 
a straightforward induction on a (using Lemma T2. 211 . 

Lemma 2.3. Let <fi be as in the proof of Theorem 1 1.21 fl) and n>2. Then, for all 
a e Z+ and b e N, we have (p( +i)(n-i)+6,o) = 0. 

We now prove Theorem II. 21 (2). 

Proof of Theorem \l./!\ (2). With Ni,N 2 as in Theorem 11.21 we suppose Ni = n, 
N2 = m < q. As before, the first differential equation for y gives us 

(2.4) y n = -a„_iy n _i aij/i - a a { e K 

while the second one for y q (after dividing through by y q ) produces the equation 

Pm,q + b m -iPm-l,q H h b = &i € K, 

Differentiating the original linear differential equation for y, k times, produces linear 
equations y n+k = L k (yi, . . . ,y n -i) in terms (over K) of 2/1, . . . , y n _i like (|2.4|) 
above. If we also differentiate the equation for y q , k times, we will arrive at another 
equation for the variables yc 

m_1 k fk\ 

P?n + k,q :=p m+ k,q + ■ ) (D J h) Pk-j+i,q = 0. 

i=0 j=0 

It is clear that the leading homogeneous forms of the P m +k,q (with respect to 
the grading above) are p m +k,q- Let <j> be as in the proof of Theorem 11.21 (1), and 
let P m +k,q denote the polynomials produced by substituting the linear forms Li 
for the variables y n +i (i — 0,1,...) into the polynomials, P m +k,q- If Pm+k,q : = 
4>(pm+k,q) 7^ 0, then the leading homogeneous form of P m +k,q is p m +k,q because we 
are substituting linear polynomials with strictly smaller degree (corresponding to 
the grading). 

Consider the following system of equations (recall that q > m and n > 2), 

(2-5) Vrn.q = 0, Pm+l,q = 0, . . . , _P(g+l) („- 1) - 1,, = 0. 

We claim that (0, ... ,0) is the only solution to (|2.5I) . Suppose, on the contrary, 
that (yi, ... , y„-i) ^ (0, . . . , 0) is a solution to l|2.5|l . and let r € {1, . . . , n — 1} be 
the largest integer such that y r ^ 0. Also, choose t G {1, . . . , q} minimial such that 



(2-6) ptr.t = 0, ptr+l,t =0, . . . , p (t+ i) r _ M = 0. 

Clearly t ^ 1, as then p r ^\ = y r = 0, a contradiction. Applying Lemma 12.21 with <fi 
(and maximality of r), examine the equation, 

, nw . /(t + l)r-l\ . 

(2-7) P(i+l)r-l,(=P(t+l)r-l,(-H H \y r ptr~l,t-l- 

Using Lemma 12.31 (with a = t — 1) and the maximality of r, we have ptr+6,*-i = 
for all b E N. Consequently, l|2.7|l and l|2.6|l imply that pt r -i,t-i = 0. Repeating 
this examination with P(t+u r —2,uP(t+i)r-3,t> ■ ■ ■ >Ptr,t ( m that order) in place of 
P(t+i)r~i,t on the left-hand side of (|2.7|l . it follows that p tr -i,t-i = for i = 1, . . . , r. 
This, of course, contradicts the minimality of t and proves the claim. 
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It now follows from Corollary II .61 that the variety determined by 

j-fm,g = 0, . . . , P(q+\){n-l)-l,q = j 

is zero-dimensional. Therefore, by an elementary dimension argument, there exists 

a subset of at most n — 1 of these equations that generates a zero-dimensional ideal. 

The degree bounds in the theorem now follow by noticing that each polynomial, 

Pm+k,q, has total degree at most m + k, so that Bezout's theorem gives us the 

bound of (<?( "- 1)+ "- 2)! 
oouna 01 ^rripTji ■ 

Finally, we remark that the methods above still hold if the characteristic of K 

is larger than (g + l)(n — 1) — 1. □ 

3. Tighter Degree Bounds 

It is possible to refine the degree bounds in our main theorem if we restrict our 
attention to fields of characteristic zero (see Remark lrt . 21 concerning the perfect field 
case) . 

Theorem 3.1. Assuming the hypothesis as in Theorem 11.21 with K having char- 
acteristic zero, the degree of the polynomial for D^y/y (J = 1, . . . , Ni — 1) over K 
in (1) is at most { N2 ^_^ 2 ) and for (2), at most ( q(Nl N \ ) ^ 1 ' 2 ) . 

Proof. We prove the result for (1) as the other case is the same. Let Ni = n, N2 = m 
and set P m +k £ K[yi, ■ ■ ■ , Un-i] {k = 0, . . . , n — 2) to be the polynomials in (|2.2() 
after substitution of the linear forms, y n +i — Li{yi, . . . , y n _i). Corresponding to 
the grading iv(yj) — j, the weight of each monomial in P m +k is less than or equal 
to m + k. Let S be the set of all (affine) solutions with coordinates in K to the 
system {P m+ k = 0}fe=o ■ Om- first goal is to bound the cardinality of S by ( m+ "7 ) ■ 

Suppose that {y^i, . . . , yi lS } is the list of all s distinct i-th coordinates of mem- 
bers of S. Since K is infinite, there exists kt S K such that yi.j ^ k% for 
j = 1, ...,s. Now, let xi, . . . , x n -i be variables and consider the new polyno- 
mials F m+ k G K[ ] produced by the substitution yt — x\ + h in the 
P m +k- As the n — 1 equations P m +k = define a zero-dimensional variety, so do 
the n — 1 equations F m+ k = 0. 

Let S denote the set of all solutions with coordinates in K to the system {F m+ k = 
0}fc=o- S mce the total degree of each F m+ k is just m + k, we have by Bezout's 
theorem, 

\s\< im ( + n -* )l = (n-iy.( m + n ; 

(m — iy. \ n — 1 

Consider the (set-theoretic) map ip : S — > S given by (xi, . . . ,x n _i) 1— > (xi 
ki, . . . , x"z\ + fcn-i)- It is easy to see that 



(3.1) EI^WH 5 ! 



Let (y\, . . . , y n -i) G S. By our choice of ki, the polynomial hi{xi) = x\ + ki — yi 
has precisely i distinct zeroes. These i roots are distinct since characteristic zero 
implies that gcd(/i,, h^) — 1. Hence, | , _1 (s)| > (n — 1)! for all s £ S, and so from 

|5|(n-l)!<|5|<( m +^- 2 )(rz-l)!, 
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we arrive at the desired bound on \S\. 

Finally, suppose that gi(x) £ K[x] is the irreducible polynomial for j/i over K. 
Since K has characteristic zero, this polynomial has distinct roots. Thus, there are 
deg(<7i) distinct embeddings a : K(yi) — > K that are the identity on K. Moreover, 
each of these homomorphisms extends to an embedding a : K — > K p. 233]. 
In particular, the deg(gi) points, (cn/i, . . . , ay n -i), are all distinct elements of S. 
Thus, from our bound on \S\, we must have 

/ \ , a, frn + n — 2\ 
deg( 5l ) < \S\ < ^ n i J. 

This completes the proof of the theorem. □ 

Remark 3.2. Examining the proof of Theorem 13. II it should be clear that the same 
result holds when K is an infinite perfect field and the hypotheses of Theorem 1 1.21 
are satisfied. 

We should also note that the proof above generalizes to bound the number of 
distinct solutions to certain systems of equations. Specifically, we have the following 
interesting fact. 

Theorem 3.3. Let w(yj) — j be the grading as above and let K be a field of 
characteristic zero. Let m 6 Z+ and suppose that {F m +k(yi, ■ ■ ■ , JM-i) = 0}^Zq is 
a zero-dimensional system of polynomial equations over K such that each monomial 
in F m +k has weight less than or equal to m + k. Then, this system will have at 
most ( m +-i 2 ) distinct solut ions with coordinates in K. 

In principle, the number of solutions for a generic system with conditions as in 
Theorem l3 . 3l can be found by a mixed volume computation and Bernstein's Theorem 
(see PP, for instance). This approach, however, seems difficult to implement. 

4. Applications to Nonlinear Differential Equations 

In the proof of Theorem 11.21 it is clear that the important attributes of the 
recursions as in (|2.1|l are that they reduce the degree and are polynomial in nature. 
In particular, it was not necessary that they were linear. For example, the system, 

yy"' + a(y') 2 +by 2 = 0, 
(l/y)" + c(l/yY + d(l/y)=0 

gives us the recurrence 7/3 +ay\ + b — 0, which has y% expressible as a polynomial in 
2/i,2/2 with strictly smaller weight. Repeated differentiation of this equation, pre- 
serves this property. In general, let h G K [zi, . . . , z n ] be a homogeneous polynomial 
(with respect to total degree) such that each monomial z a — z" 1 ■ ■ ■ z"™ has 

n 

(i — l)ct; < n. 

i=l 

If the hypothesis of Theorem 1 1.21 are weakened to allow y to satisfy an equation of 
the form, D n y = h(y, Dy, . . . , D^^'y), then the proof applies without change. A 
generalization along these lines was also considered by Sperber in |Jj , however, the 
techniques developed here give us degree bounds just as in Theorem ll.2l and also 
allow for a positive characteristic. 
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